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Abstract 

The pseudoscalar resonance or “A-funnel” in the Minimal Supersymmetric Standard 
Model (MSSM) is a widely studied framework for explaining dark matter that can yield inter¬ 
esting indirect detection and collider signals. The well-known Galactic Center excess (GCE) at 
GeV energies in the gamma ray spectrum, consistent with annihilation of a < 40 GeV dark matter 
particle, has more recently been shown to be compatible with significantly heavier masses following 
reanalysis of the background. In this paper, we explore the LHC and direct detection implications 
of interpreting the GCE in this extended mass window within the MSSM A-funnel framework. We 
find that compatibility with relic density, signal strength, collider constraints, and Higgs data can 
be simultaneously achieved with appropriate parameter choices. The compatible regions give very 
sharp predictions of 200 — 600 GeV CP-odd/even Higgs bosons at low tan/3 at the LHC and spin- 
independent cross sections « 10“^^ pb at direct detection experiments. Regardless of consistency 
with the GCE, this study serves as a useful template of the strong correlations between indirect, 
direct, and LHC signatures of the MSSM A-funnel region. 
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I. INTRODUCTION AND MOTIVATION 


The Galactic Center (GC) of the Milky Way galaxy is the densest dark matter region 
in our vicinity and has long been earmarked as the most promising target for searches of 
dark matter (DM) signals. Intriguingly, recent years have seen a persistent and statistically 
signihcant excess in the gamma ray spectrum peaking at 2 —5 GeV originating from the GC, 
above what is predicted from known sources and conventional astrophysics [M2]. The signal 
was initially reported to be compatible with ~ 40 (10) GeV dark matter annihilating into bb 
(rr), with an annihilation cross section (av) ~ 0(10“^®) cm^/s. Since this is approximately 
the annihilation cross section expected of a thermal relic, a dark matter interpretation of 
this excess presents itself as a very tantalizing possibility. This prospect has been explored 
by many authors in various contexts (see, for instance Refs. [91201 El E] and references 
therein), including the Minimal Supersymmetric Standard Model (MSSM) [TSl - ITS] . More 
recently, it has been shown that this excess might be attributable to unresolved point sources 
[TUH?T] . although a conclusive verdict has not been reached. 

Recently, the Fermi-LAT Gollaboration has presented an analysis of the region around 
the GG with four different variants of foreground/background models, hnding, for every 
variant, significant improvements in the agreement with data when an additional component 
centered at the GG with a peaked profile (NFW, NFW-contracted), i.e. a dark matter-like 
spectrum, was included in the fits [ 221123 ] (see also Ref. [ 12 ] for an attempt at accounting 
for systematic uncertainties in the background). From a dark matter perspective, a recent 
study na found these additional components for the four choices of background models to 
be compatible with several annihilation channels {WW, ZZ, hh,tt) and significantly higher 
DM masses (165 GeV for bb, 310 GeV for tF) than previously thought possible. Similar 
conclusions were also reached in Refs. ini and [H], which reported that a higher mass 
(175 — 200 GeV) dark matter annihilating into ti could give reasonable hts to the signal. 

This relaxation of the allowed range of dark matter masses compatible with the GG ex¬ 
cess (GGE) has particularly interesting implications for MSSM dark matter, as it opens 
up the possibility of explaining the signal with the well-known pseudoscalar resonance or 
“A-funnel” mechanism, where the dark matter relic density is set by resonant s-channel an¬ 
nihilation through the pseudoscalar A, with rriA ^ 2m^ (y represents the lightest neutralino, 
which is the dark matter candidate). The pseudoscalar resonance has been studied in con- 
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nection with the GCE outside the MSSM in Refs. [23H26]; however, realizing the mechanism 
in the MSSM is of particular interest given that the MSSM remains one of the most familiar 
and widely studied Beyond the Standard Model (BSM) theories. Previous hts to the GCE 
with < 50 GeV did not allow for this possibility in the MSSM due to constraints on niA 
from direct LHC searches |23 [28] (although this constraint can be circumvented in the the 
Next-to-Minimal Supersymmetric Model (NMSSM), allowing for an NMSSM explanation of 
the GCE [25],|29]). This incompatibility is lifted if, as discussed in Ref. |T5], < 165 (310) 

GeV annihilates into bb {tt), allowing for large enough to evade collider constraints. 

The aim of this paper is to explore whether, given this wider range of allowed masses, 
the MSSM pseudoscalar resonance can give reasonable hts to the GCE, consistent with 
stringent constraints from relic density, indirect/direct detection, collider search limits, and 
Higgs data. Since the mechanism requires a light (~ 200 — 500 GeV) pseudoscalar, the 
SM-like nature of the 125 GeV Higgs boson is particularly constraining as the heavier CP- 
even Higgs is at the same mass as the pseudoscalar and can mix with the 125 GeV Higgs, 
resulting in deviations from SM-like properties inconsistent with measurements. For such 
light, non-decoupled heavier Higgs bosons, the Higgs sector needs to be “aligned” [301 - 135] to 
maintain SM-like properties for the 125 GeV mass eigenstate. As we will show in this paper, 
this can indeed be achieved while simultaneously satisfying all other DM requirements. 

A successful realization of neutralino dark matter along with the GCE through the pseu¬ 
doscalar resonance requires very precise choices of parameters in order to simultaneously 
achieve resonant annihilation, the Higgs mass, and alignment in the Higgs sector (this is 
also the reason why extensive scans in the MSSM parameter space pMH] fail to uncover it 
as a viable explanation of the GCE). It is nevertheless worthwhile to pursue this direction 
for several reasons. First, the A-funnel is one of several “traditional” mechanisms in the 
MSSM that have been widely studied for a long time, and its compatibility with a possible 
DM signal is therefore of considerable interest. Second, while most scenarios put forward 
to explain the GCE could potentially be constrained by stringent spin-independent direct 
detection limits (indeed, avoiding these limits itself involves some nontrivial hne-tuning of 
parameters in supersymmetric models [361138] L the A-funnel naturally gives small direct 
detection cross sections and is automatically safe from these bounds. Most importantly, the 
framework is eminently predictive, giving very specihc predictions for heavy Higgs bosons 
that will be probed at the 13 TeV LHC and future colliders, as well as direct detection 
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cross sections that may be probed by the next generation of experiments. Independent of 
these considerations, and independent of the applicability to the GCE, this study serves as 
a valuable template of the conditions necessary for the existence of a light pseudoscalar in 
the MSSM together with indirect detection signals of dark matter via the A-funnel. 

The outline of the paper is as follows. Section II introduces the parameter space relevant 
for the study and discusses dark matter aspects such as the annihilation cross section and 
relic density. Section III is devoted to a discussion of various constraints from direct detec¬ 
tion, indirect detection, collider constraints, Higgs data, and vacuum metastability. Section 
IV presents the details of our scans and the best £t regions to the GCE. Predictions for 
the 13 TeV LHC and future direct detection searches are presented in Section V. We sum¬ 
marize our results in Section VI. The Appendices contains additional details on the MSSM 
parameters and hts to the GCE. 

II. THE MSSM PSEUDOSCALAR RESONANCE: DARK MATTER ASPECTS 

In i?-parity conserving supersymmetric models, the lightest supersymmetric parti¬ 
cle (LSP) is stable. If it is also neutral, it can be a dark matter candidate. In the MSSM, 
the LSP is often assumed to be the lightest of the neutralinos, the neutral superpartners 
of the gauge bosons and Higgs bosons (Bino, Wino and Higgsinos respectively). The Wino 
and the Higgsinos tend to annihilate too efficiently to explain the observed dark matter 
abundance. However, the Bino can yield the correct relic density via various mechanisms, 
including resonant annihilation via the pseudoscalar, and has long been regarded as the 
favored dark matter candidate. 

We perform our study in the phenomenological MSSM (pMSSM) [35], which is dehned 
in terms of 19 parameters, which are taken to be independent at the weak scale. Of these, 
our analysis will be entirely determined by the following seven parameters: 

• Ml, the Bino mass parameter. The dark matter is mostly Bino, so this is also approx¬ 
imately the mass of the dark matter candidate ~ Mi. 

• fi parameter. This is the Higgsino mass, and controls the Higgsino fraction in the dark 
matter particle y. As we will see later, the relic density, signal strength, and direct 
detection cross section all depend sensitively on this fraction. 
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• tan/?, the ratio of the up- and down-type Higgs vacuum expectation values (vevs). 

• rriA, the heavy Higgs mass. This is the mass of the pseudoscalar that mediates the 
resonance (hence niA ~ 2m^) as well as the mass of the heavier scalar, which feeds 
into Higgs phenomenology and expected direct detection cross-sections. 


mQ^,mu 3 , right handed stop masses, which contribute signihcantly to the 

mass of the observed 125 GeV Higgs boson. In this paper we take the stop mass scale 


Ms = = rnls- 


• At, stop trilinear coupling. This determines the mixing in the stop sector and is again 
a relevant parameter for the mass of the observed Higgs boson. 


All other masses, such as the other gaugino (wino and gluino) and sfermion masses, are 
assumed to be heavy and decoupled from the analysis. 


A. Dark Matter Composition 


The lightest neutralino in the MSSM is a combination of the Bino, Wino, and neutral 
Higgsinos: 

X = NnB + NuW + . (1) 

As mentioned above, we are mainly interested in the region of parameters where the lightest 
neutralino is predominantly a Bino, hence An ~ 1, A 12 = 0, and Ai 3 ,Ai 4 1. In this 
regime, the Bino mass parameter Mi and the neutralino components are approximately | 


Ml = rriy + 
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B. Relic Density and Signal Strength 


Both the relic density and the present day annihilation cross section are driven by the 
process yy —)■ // with the psendoscalar A in the s-channel (we are interested in the case 
where the fermion / is either bort for compatibility with the GCE). When the process occnrs 
close to resonance, it is well-known that the annihilation cross-section in the early nniverse 
(which sets the relic density at the time of freeze-out) is substantially different from that at 
present times (which sets the signal strength htting the GCE) due to thermal broadening of 
the resonance during the former stage |1D]. Thus, with appropriate parameter choices, one 
can scale the relic density and the present annihilation cross section independent of each 
other, thereby achieving better agreement with both measurements; this degree of freedom 
is not afforded in non-resonant scenarios, where these two quantities are strictly related to 
each other. 

To understand this interplay, consider a simplihed model describing a Majorana DM 
particle y coupled to a psendoscalar A through the interaction Lagrangian 

-Dint = Waxx^X'l^X + WaffAf'y^f. (3) 


The entire parameter space of the model is then determined by mA,my^,yaxx i/aff- A 
crucial parameter in our analysis is the degeneracy parameter 


(5 = |1 - 


(4) 


which characterizes the proximity to the resonant regime. We are interested in scenarios 
where 5 ~ 0. 

The resonant annihilation cross-section at a given temperature T is 


(o'v) ~ ^ Q ^ 

y/nm\TA 

where x = m^/T and T^ is the decay width of A, 
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where x/ is the value of x at freeze-out. This expression can be rewritten in a more illumi¬ 
nating form as [2S] 
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Likewise, the DM annihilation cross-section today is 
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Assuming that mA ~ 2m^ so that the second term dominates in the denominator, one 
obtains (for 2 m^ < nia) 
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Gomparing Eq.j^and Eq. IT, it is clear that the relic density and the current annihilation 
cross-section can be independently scaled with judicious choices of Uaff and yaxx\J 
terms of the fundamental MSSM parameters, these couplings are given by: 
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( 12 ) 


where v = 174 GeV and gi is the SM U{1)y gauge coupling. Note from the above that a 
non-vanishing Haxx coupling requires a non-vanishing Higgsino component in y. From the 
expressions for Vn, V 13 , listed previously, we thus see that, for given values of ttia and 
tan/?, the desired relic density and an annihilation cross-section consistent with the GGE 
can be obtained simultaneously by appropriately choosing /i and 5 (equivalently, m^). 


III. CONSTRAINTS 

As mentioned in Sec. the relevant A-funnel parameter space is constrained from several 
directions. Higgs phenomenology in our set-up is very directly linked to the GGE, hence 
LHG direct searches as well as the properties of the observed 125 GeV Higgs put stringent 
constraints on this scenario. Gonsistency with all collider observables can then create tension 
with constraints from requiring the stability of the electroweak vacuum. In addition, since 
the GP-even heavy Higgs H is expected to be approximately degenerate in mass with A, 
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contributions to the spin-independent direct detection cross-section from iif-exchange might 
be relevant. Finally, there are also several current and future indirect detection experiments 
that can probe the process of interest in this paper. In this section we detail the current 
status and future prospects in all of these different directions. 

A. Collider and Higgs Sector Constraints 

In the absence of CP-violation (which we assume in this paper), the physical spectrum 
of the Higgs sector consists of two CP-even Higgs bosons, h and iP, one CP-odd state A, 
and a pair of charged Higgs bosons, H^. Direct searches for these heavier Higgs bosons at 
the LHC rule out a signihcant part of parameter space. ATLAS and CMS direct searches 
for charged Higgs bosons [2ZIEH] rule out mH+ < 160 GeV (recall that = m\ + rri^ 
at tree level). Likewise, there exist strong limits from searches for A/H —)■ rr|H], which 
provide the strongest limits, although these depend on tan (3 and can be evaded for small 
values of tan/3.^ 

Beyond these direct constraints, a small ttia is still in tension with Higgs data, as a light 
CP-even Higgs (m^ ~ rriA in the MSSM) tends to mix with the 125 GeV state and cause 
deviations from SM-like properties. This is a particularly strong constraint in our framework 
and dictates what values our parameters can take, hence we will now study this constraint 
in some detail. 

The MSSM Higgs sector consists of two doublets, Hu and Hd] the former couples to all 
the up-type fermions and the latter to the down-type fermions and charged leptons. The 
neutral components acquire vacuum expectation values Vu and Vd with tan/3 = tp = Vu/vd 
and = 174 GeV. One can dehne a “Higgs-basis”, where a single held acquires 

all the vev: 


Hsm — spHu + cpHd, (13) 


^ Light niA/mH and heavily mixed stops (as usually needed for a 125 GeV Higgs in the MSSM) can also 
give large contributions to various flavor observables, for example Bg and B -A- Xgj. However, 

in this work we will mainly be interested in moderate to small value of tan (3, hence there is no large 
enhancement of these effects. Moreover, the size of these contributions are heavily dependent on the signs 
of various contributions (see e.g. Ref. [HI), and consistency with all measured values could be obtained 
by tuning such cancellations. 




Hnsm — —cpHu + SjsHd, ( 14 ) 

where = sin/3, C /3 = cos/3, (Hsm) = v, and (Hnsm) = 0. The couplings of these states 
to the SM fields are: 

dd/uu/VV SM 

Qr = Q 

jHsm ^ 

9Zs.. = -9^^/t^ 9 Zsm = 0 > ( 15 ) 

where VV, uu, dd refer to all vector, up-type and down-type states respectively, and 
refers to the SM value of these couplings. Note that there is no coupling between the 
states and the or between the gauge bosons and 
The mass eigenstates, h and H, can be written as mixtures of the Higgs basis helds, 

H = ZmHsm + k,nsmHnsm , ( 16 ) 

where Zsm — ~>^sm — = cos(q: — /3) and = h^sm = ^a-/3 = sin(Q; — /3), and a is 

the angle of rotation from the (hf„, Hd) basis to the mass eigenstates. We want to identify 
the lightest CP-even mass eigenstate, h, with the recently observed 125 GeV scalar; given 
that all measurements suggest that its properties are SM-like, we also want to identify it as 
the SM-like field in the Higgs basis. That is, we require 

hi25 = h K, Hsm- ( 17 ) 


This requirement of vanishing mixing between the Hnsm state and the 125 GeV Higgs, 
corresponding to Zsm ~ 0; can be rewritten in terms of the fundamental parameters as 

[301 EH 
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(18) 


where Ms is the geometric mean of the stop masses and 


Xt = At — , Yt = At + gtfd. (19) 

Note that when the second Higgs becomes heavy {mH » frih), this relation is automatically 
satisfied; this is the familiar decoupling effect. Otherwise, one requires alignment without 
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decoupling [30l[3T], brought about by an accidental cancellation in the fundamental param¬ 


eters of the theory so as to satisfy Eq. 18 For small and Ms ~ 0(1) TeV, large values 


of At/Ms are required to obtained an experimentally consistent Higgs mass whereas large 


values of {^At)/Ms lead to close to alignment conditions [3011^. 


The CMS and ATLAS collaborations present both the precision measurements of the 125 
GeV Higgs and the searches for H —)■ WW/ZZ as ratios to the expectations from a SM 
Higgs of the same mass. The predicted rate at the LHC for the decay of the mass eigenstate 
i = {h, H} into some hnal state XX as a ratio to the SM value is given by 


Kix = W'/SM) X {BFf:,^/SM) . 


( 20 ) 


where SM in the denominators denote the corresponding values for a SM-like Higgs of the 
same mass. For a 125 GeV SM-like Higgs, the dominant decay mode is into a pair of b- 
quarks (~60%), followed by WW; hence the total width is dominated by the width into b 
quarks. The largest deviation from mixing effects is expected in the precision measurements 
of h —)■ WW. This number is reported to be 'R-ww ~ by ATLAS [13] and R-ww ~ 

0.83 ± 0.21 by CMS [Hj. In our analysis we will take a conservative approach of assuming 
that observational consistency is obtained (that is, the Higgs sector is sufficiently aligned) 
for R-ww between 0.7 —1.3. This range will narrow with additional data, and measurements 
at the level of 10% are expected at the high luminosity LHC |35l H6] . 

B. Vacuum Metastability 

Another important constraint on these parameters comes from vacuum metastability. 
Large values of the soft stop trilinear coupling At, required for the Higgs mass and alignment 
(discussion above), can result in the appearance of charge- and color-breaking minima in 
the scalar potential of the MSSM. The condition for either these minima to be energetically 
unfavorable or the tunneling to these minima to have lifetimes longer than the age of the 
Universe leads to the approximate bound m 



( 21 ) 


where = uiq^ + + /i^, and r = In our analysis we assume 

rugg = m ^3 = M|, so that r = 1. Minimization conditions of the Higgs potential give 
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ml = m\cos^ (3 + 0.5m| cos(2/3), hence the condition for vacunm metastability can be 
written as 

Af < 6.8M| + 60m^cos^/5 + 30 m| cos(2/d). (22) 

It is worth keeping in mind that this is only an approximate bound and depends on several 
assumptions (see Ref. m for details). However, consistency with the above provides a rough 
guide for the feasibility of the parameter region under investigation. 


C. Direct Detection 


Direct detection possibilities focusing on the H-funnel in the MSSM have been studied 
in Refs. |l8ti50] . The pseudoscalar A does not mediate spin-independent WIMP-nucleon 
scattering. Instead this cross section asi comes from light and heavy CP-even Higgs boson 
exchanges in the t-channel, facilitated by the Bino-Higgsino mixture of the LSP necessary to 
obtain the correct relic density. There are also contributions from tree level squark exchange 
in the s-channel and from gluon loops [HU |52] , but these are negligible when the sfermions 
are heavy. The cross section then depends only on Mi, tan/3 and /i. 

For given values of mA and tan /3, requiring the correct relic density and GCE leaves no 
free parameters, thereby hxing the direct detection cross section. This cross section in our 
region of interest can be written as approximately ^5] 
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(23) 


where Fu ~ 0.15, ~ 0.13 (the up and down type quark content respectively of the 

nucleon, proton or neutron), tjs = tan/3, m^ is the mass of the nucleon, and 
is the reduced mass. For the correct dark matter relic density obtained via the H-funnel, 
this cross section is generally around 10“^^ pb [1811501153] . well below existing bounds from 
XENONIOO [5l] and LUX [55], which currently rule out asi ^ 5 x 10“^° pb. Note that while 
the annihilation processes that determine the relic density as well as indirect detection signals 
are s-channel and therefore enhanced by the resonance, the direct detection cross-section is 
mediated by t-channel processes and does not receive this enhancement. Such small direct 
detection cross sections are therefore a generic feature of this region of parameter space. 
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Crucially, this cross section still lies above the neutrino background and is therefore within 
reach of future detectors, although detection will still be challenging. 

As is well-known, an exception to this generic feature can occur for negative values of the 
p parameter due to destructive interference between the light and heavy Higgs exchange con¬ 
tributions, giving cross sections several orders of magnitude below the neutrino background 
cross section |l8l |50]. Such blind spots can in general occur at any dark matter mass, but 
their appearance in the A-funnel framework is more strongly constrained as we also need 
niH ~ ttia ~ 2m^. Approximating the up- and down-type quark content in the nucleus 
as roughly equal, this cancellation condition in the A-funnel region can be formulated as 
approximately |1H] 

rriA ~ tan/?^ . (24) 


With TeV scale values of p necessitated by relic density constraints and 0{1) values of tan (3 
required by collider constraints (see Sec. Ill A), Eq.[^ implies that the cancellation can only 
occur for large ^ 650 GeV, beyond the mass range of interest from the point of view 
of the GeV excess. Hence all parameter combinations of interest should predict a small but 
tractable (~ 10“^^ pb) direct detection cross section (we will see in the subsequent sections 
that this in indeed realized, see Fig.[^. 


D. Indirect Detection 

Gurrently the strongest bounds on the annihilation cross section are given by the 
Fermi/LAT analysis of 6 years of data on 15 known dwarf galaxies [56]. For 100 — 300 
GeV dark matter, which is our region of interest, this analysis constrains the annihilation 
cross-section to be less than ~ a fewxlO”^® cm^/s. The cross section required to explain 
the GGF is also in this region over this mass range (see US]), hence the dwarf constraints 
are in some tension with a DM interpretation of the GGF. However, the large uncertainties 
in the dark matter distribution (J-factor) in these dwarf galaxies leave room for compatibil¬ 
ity (see Fig. 8 in Ref. [56]). For instance, the 95%G.L. annihilation cross-section exclusion 
limit for a 100 GeV WIMP annihilating to bb is 2.2 x 10“^®cm^/s and has a Icr error inter¬ 
val of [9.0 X 10“^^, 5.6 X 10“^®] cm^/s, which is compatible with the cross section interval 
[3.1 X 10“^^, 8.8 X 10“^®] cm^ needed to £t to the GGF at this mass. A signal was reportedly 
seen in the new dwarf galaxy candidate Reticulum H [57], found in the first year DFS data 
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[58] . consistent with a dark matter of mass ~ 40 — 200 GeV annihilating into bb with a 
cross section (av) ~ 10“^® cm^/s, although this was later found to be inconsistent with the 
new PASS 8 diffuse emission model used to analyze Reticulum II [SS]- Bounds similar to 
those from the Fermi dwarf observations are also found by the Planck satellite from CMB 
measurements [60] , 

Likewise, since DM of interest in this paper annihilates primarily through hadronic chan¬ 
nels {bb and tF), this is expected to generate a signihcant flux of antiprotons. There already 
exists some tension between models that explain the GCE and derived constraints from an¬ 
tiproton bounds on dark matter annihilation [MHM]- However, calculation of the antiproton 
flux suffers from signihcant uncertainties related to the propagation model in the galaxy (see 
[ 63 H 66 ] and references therein), and the GGE can be made compatible with the measured 
antiproton hux for conservative choices of propagation model parameters. 

Bounds on the dark matter annihilating cross-section into quarks are also obtained by 
neutrino experiments like IceGube. The most current results from the IceGube-79 experi¬ 
ment exclude {av) > 2 x 10“^^ cm^/s into bb at 90% conhdence level [SZ]. This lower limit 
is ~ 10^ larger than the cross-section required for the GGE [15] and thus irrelevant. 

Therefore, no indirect detection results robustly rule out a DM interpretation of the 
GGE at present, although future measurements, particularly from Fermi-LAT observation 
of dwarfs, AMS-02 antiproton results, and the GMB could have interesting implications. 


IV. NUMERICAL RESULTS 

Building on the parameter space and constraints described in the previous sec¬ 
tions, we present the hts to the GGE excess in this section. We used the following 
tools for our numerical analysis: the neutralino relic abundance and annihilation cross- 
section was calculated with Micromegas-4.1.7 [6H], the MSSM particle spectra were 
computed using SuSpect-2.41 [69], and the Higgs phenomenology was obtained with 
FeynHiggs-2.11.0 [701 - I74] . 

For the gamma ray spectrum corresponding to the signal, we follow the approach em¬ 
ployed in Ref. [IS] and consider two of the four spectra presented in Fig. 13 of Ref. [23]^. 

^ The first version of our paper used the spectra presented in Ref. |^, and Ref. [23j is the corresponding 
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which were derived by fitting the excess over various choices of background as exponentially 
cut off power laws (see Ref. [I51E2] for further details). The four spectra are referred to 
as spectra (a)-(d) in Ref. [U], and just as they do, we pick spectra (b) and (d) for our 
analysis; spectrum (a) is very similar to what has been studied for light (m^ < 40 GeV) DM 
in previous papers and not amenable to the MSSM, whereas spectrum (c) is very similar to 
spectrum (d) and does not yield any new insight. 

Spectrum (b) corresponds to a fit with OB stars as cosmic ray (CR) sources and a tuned 
index for pion production within the solar circle (see [221123]); the analysis in Ref. [T^ found 
it to be well fit by 75 — 95 GeV DM annihilating into bb or < 200 GeV DM annihilating 
into ti. Annihilation into gauge or Higgs bosons were also found to give good fits, but these 
are irrelevant for our analysis since they are always subdominant channels in the MSSM 
pseudoscalar resonance scenario. Note that spectrum (b) is also in agreement with other 
studies performed in Refs. ini and [18], which also found that 175—200 GeV DM annihilating 
into ti could be compatible with the GGE. Likewise, spectrum (d) corresponds to a fit with 
OB stars as cosmic ray (GR) sources but with only the intensity of pion production tuned 
(using pulsars instead of OB stars gives a very similar spectrum); Ref. [I5] found it to 
correspond to higher mass DM, with 130 — 165 GeV DM annihilating into bb or 250 — 310 
GeV DM annihilating into ti giving good fits. 

In this section, we will perform fits to the two spectra (b) and (d) with the idea of gaining 
intuition about the range of possibilities that the GGE allows for the MSSM pseudoscalar 
resonance. We note that the continuous region spanning spectra (b) and (d) could also 
plausibly explain the GGE for some reasonable background, but do not pursue this direction 
any further. 

A. Fit Procedure 

The astrophysical information regarding the distribution of dark matter is encoded in the 
J-factor 



(25) 


publication that recently appeared; we have chosen the spectra from Ref. [23] that correspond most closely 
to the spectra we used in the first version. 
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where Al] is the region of interest (ROI), l.o.s. stands for line of sight, and p is the dark 
matter density. Jean. = 2.0 x lO^^GeV^/cm® is the canonical value of the J-factor obtained 
from evaluating the integral with an NFW prohle. Following the analyses in Ref. [15], we 
parametrize the uncertainty in the dark matter density prohle with the factor which is 
allowed to vary between [0.14,4]. 

The gamma-ray spectrum is computed for the following MSSM parameters: 

• The pseudoscalar mass is allowed to vary over 200 GeV < mA < 700 GeV. Below 

200 GeV, we hnd that the Higgs sector cannot be sufficiently aligned while remaining 
consistent with bounds from H/A ^ from the 8 TeV LHG run. We terminate 

the scan at 700 GeV since good hts to the GGE (either spectrum (b) or (d)) are not 
expected for > 310 GeV. 

• tan /3 is scanned over the range 4 < tan (3 < 10. Below tan /? = 4, extremely heavy 
(multi-TeV) stop masses are required to reproduce the Higgs mass, and large log 
resummations become important. Above tan/3 ~ 10, uia ^ 350 GeV is inconsistent 
with the LHG H/A ^ bound. Masses heavier than this do not give good hts to 
the GGE since rriA ^ 310 GeV (m^ > 165 GeV) requires annihilation primarily into 
tt, but for large values of tan f3 the leading annihilation channel for the pseudoscalar 
is into bb. 

• For given values of and tan /3, we next scan over 6 (equivalently, as shown in 
Eq. 1^ and p for points such that 

— the relic density constraint is satished: the neutralino makes up all of dark matter 

(0.08 < < 0.16); and 

— the annihilation cross section (av) is within the 2a best-ht annihilation cross- 
section contours from Ref. [T5] . 

We scan over S G [0, 0.1] in order to stay close to resonance, and over p G [0.7,10] TeV 
in order to obtain a mostly bino DM. 

• Next, we scan over the stop masses {Ms = mg^ = m^g) G [0.7,12.7] TeV and the stop 
trilinear coupling G [5, 25] TeV for points satisfying 

— 122 < ruh < 128 GeV; and 
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alignment in the Higgs sector. 


We take the branching ratio to WW normalized to the SM value 'R-ww be a measure 
of alignment and select (for each niA, tan 13, /i, combination) the combination of Ms 
and At that gives R-ww closest to 1 while maintaining 122 < rrih < 128 GeV. 


• All other MSSM input parameters (gaugino/wino masses, trilinear couplings, slep- 
ton/squark masses) are set to 5 TeV so that they decouple from this analysis. 


The goodness of £t is obtained by performing a analysis between the gamma-ray 
spectrum obtained from Micromegas and the GCE (Fermi spectra (b) and (d)). For a given 
MSSM point, the is calculated as: 


=i: 




at 


(26) 


where the subscript k runs over the 20 energy bins of the Fermi/LAT measurement 
dN/dE is the gamma-ray spectrum obtained from Micromegas, the subscript obs denotes 
the spectrum consistent with the Fermi excess (he. spectrum (b) or (d)), cxfc denotes the sta¬ 
tistical uncertainty lEI. and ^ is the value of Jh G [0.14,4] that minimizes the value. The 
analysis includes statistical errors, but neglects possible systematic errors from modeling 
backgrounds near the Galactic Genter. 


B. Fit Results 

The hts resulting from the above procedure are presented in Fig. [T]as contours of x^ in 
the m^-tan/3 plane for Fermi spectrum (b) and (d). The pink crosses in each panel denote 
the points with the best fit to the corresponding spectrum; the gamma-ray spectra of these 
best £t points are presented in Fig. along with the MSSM parameters ^ In Fig. [T]we also 
include, in solid black lines, the 1-a and 2-a bounds from A/H ^ searches at the 8 TeV 
LHG [H]; points that lie above these curves in the shaded region are inconsistent with these 
bounds. These rr searches, however, lose sensitivity at low tan [3, hence light pseudoscalars 
can mediate DM annihilations capable of explaining the GGF in this region. The dashed blue 

^ It is worth keeping in mind that the absolute value of does not have a proper statistical significance 
without a full analysis of all uncertainties in the signal and theory prediction. 
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FIG. 1: Contours of in the m^-tan/? plane from fitting the gamma-ray spectrum from the 
MSSM pseudoscalar resonance to Fermi spectrum (b) (left panel) and spectrum (d) (right panel), 
corresponding to “OB stars index scaled” and “OB stars intensity scaled” spectra from Fig. 13 
of Ref. [23] (see Ref. [I5l[22[ |23] for further details). Red (blue) contour regions denote the best 
(worst) fits. The contours are plotted in intervals of 10 in the range 20 < < 110. The pink 

crosses denote the lowest x^ value in the scan and hence represent the best fit points. Solid black 
lines mark the l-cr and 2-fT exclusion limits (shaded region above the solid black lines excluded) 
from the negative search results for H/A —)• at the 8 TeV LHC run. Dashed blue lines denote 

contours of the ratio current Higgs data from the 8 TeV LHC favors 0.7 < ^1-3 (see 

text for details). 


lines correspond to contours of 'TZww defined in Eq. 20 TZ^yy = 1 represents a completely 
SM-like Higgs, and any mixing with the non-SM Higgs causes deviations. Current Higgs 


data from the LHC allow for 0.7 < TZ^yy < 1.3, as discussed in Section HI A This leads to 
the requirement of large p and hence small couplings [c.f. Eqs. |^and 12 between A and x- 
This generically requires close to resonance conditions 2m^ ~ rriA for consistency with both 
the GCE and relic density. 

We found that the x^ value did not change signihcantly between distinct values of (p, 6, At, 
and Ms) for the same rriA, tan/?. This is expected, since the fit quality is driven by the 
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FIG. 2: Gamma-ray spectra for the best-fit points corresponding to the pink crosses in Fig. 
The gamma-ray spectra from Mi cromegas (blue line) for the best fit points are superimposed on 
Fermi spectrum (b) and spectrum (d) (black points) on the left and right panels respectively. 
The gray band denotes statistical uncertainties (from [23]). Numerical values of the corresponding 
MSSM parameters and the leading DM annihilation channels are also listed. The value of the higgs 
mass, relic density, annihilation cross-section and spin-independent scattering cross-section {nih, 
D/i^, {(Jv), asip) for the best fit point of spectrum-b/d are (126 GeV,0.082, 3.849 x 10“^® cm^/s, 
1.689 X 10-^2 pb)/(127 GeV,0.11, 3.56 x lO-^® cmVs, 4.392 x lO-^^ pb) 


shape of the spectrum, which is controlled mainly by tan (3 via the branching ratios, and the 
position of the peak, which is controlled by mA{^ 2m^). Although the £t should also depend 
on the signal strength, which is controlled by p and S via the annihilation cross section and 
relic density, the freedom in choosing G [0.14,4], which essentially rescales the signal 
strength, smears out this dependence. In our region of interest, we find that 6 < 0.04 while 
Ms, At, and /i all take multi-TeV values; we present contour plots of these parameters in 
Fig. [^in Appendix [Aj The condition for vacuum metastability, Eq. 22, is also found to be 
satisfied in most parts of the parameter space allowed by the 8 TeV LHC A/H ^ r+r“ 
bounds (see Fig. [^in Appendix [A]) . 

From the left panel of Fig. the best fit regions to Fermi spectrum (b) appear to 
be separated into two distinct islands. The < 250 GeV region has relatively low 
for all values of tan (3. In this region, annihilation into top quark pairs is kinematically 
forbidden, so the dominant annihilation channels is bb for all values of tan (3. Recall that an 
approximately 100 GeV DM particle annihilating into bb can £t the GGE na; this region 
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reflects this behavior. However, we see that this region is incompatible with the 8 TeV 
LHC AjH —)■ r+r“ bounds and/or the Higgs data (that is, Tlww ~ region, 

signaling that the heavier CP-even scalar is so light that alignment does not work well). A 
second island opens up at 350 GeV < ^ 450 GeV, when annihilation into ti becomes 

kinematically feasible, and tan f3 < 6. This is consistent with Ref. na hnding a 200 GeV 
DM annihilating into ti providing a good £t to spectrum (b). Note that the best £t point 
occurs at the lowest allowed value of tan/3(=4) in our scan, where the coupling of A to top 
quarks is the largest. The fit deteriorates as tan/? gets larger, as the branching ratio into bb 
gets larger due to the tan 13 enhancement of the Abb coupling. This region is also compatible 
with Higgs data as 'R-ww ~ from the current A/H ^ bounds. Beyond 

this island, the ht deteriorates rapidly as ttia and/or tan/3 are increased. 

Similar patterns are observed for the £t to spectrum (d). A small region of good £t 
exists at ~ 300 GeV and low tan/3, safe from the A/H —)■ r+r“ bounds and borderline 
compatible with Higgs data. Again, DM in this region annihilates dominantly to bb since 
ti is kinematically forbidden, and this observation is compatible with Ref. [15], where DM 
with mass 130 — 165 GeV annihilating into bb was found to give good fits to the spectrum. 
A second region with better hts is again observed for larger tua once decay into ti opens 
up. This regions roughly spans 450 GeV < uia ^ 600 GeV and tan /3 < 8, and appears to 
correspond to the 250 — 310 GeV DM annihilating into ti region reported in Ref. [15] as a 
good fit to spectrum (d). Similarly to spectrum (b), the best fit occurs for small values of 
tan/3 : tan/3 ~ 4.0. This suggests that a DM candidate that annihilates signihcantly into ti 
with BR{xx ti) = 0.66 at the best ht point) provides the best ht to spectrum (d). This 
can be confirmed by comparing the shape of the spectrum in Fig.[^ right panel, which hts 
the shape of Fermi spectrum (d) quite well. Finally, the ht deteriorates for larger ttia and 
tan/3 values and we do not expect any good hts beyond the region shown in the plot. 

1. Fit to a Modified Spectrum 

So far, we performed hts to spectra (b) and (d) as dehned in Ref. [15], corresponding 
to the “OB stars index scaled” and “OB stars intensity scaled” spectra from Fig. 13 of 
Ref. [23], which were obtained by modeling the excess with an NFW profile with a single 
power law with an exponential cutoh. This mimics what is expected of a dark matter source. 
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FIG. 3: Contours of in th® m^-tan /3 plane from fitting the gamma-ray spectrum from the 
MSSM pseudoscalar resonance to the “pulsars index scaled” spectrum from Fig. 18 of Ref. [23] . 
Red (blue) contour regions denote the best (worst) fits. Black and blue contours are as in Fig.[^ 

and serves the purpose of demonstrating how the preferred theory parameter space changes 
for two different choices of interstellar emission models of the background (matching the 
philosophy in Ref. na). However, Ref. [23] also hnds signihcantly better hts to the excess 
if more freedom is allowed in the £t - in particular, if the spectrum of the NFW prohle is 
modeled with a power-law that is allowed to vary per energy band over the 1 - 100 GeV 
range; the resulting spectra for various choices of interstellar emission models are presented 
in Fig. 18 of Ref. [23] • In order to study how the MSSM £t is affected if the latter is 
used, we performed a similar £t (as described above) to the “pulsars index-scaled” spectrum 
from Fig. 18 of Ref. |23|; the result is shown in Fig.[^ We hnd that the overall fit quality 
worsens due to the tail of the spectrum, but the best £t regions in the MSSM parameter 
space still closely match those from the fit for spectrum (b) (see Fig.[^ (left)); consequently, 
the theoretical implications from htting to spectrum (b) (discussed below) will also apply 
in this case. 


20 





V. PREDICTIONS FOR THE LHC AND DIRECT DETECTION EXPERIMENTS 


A. LHC Prospects 

There are several projections for the 14 TeV LHC provided by the CMS and ATLAS 
collaborations for heavy Higgs searches in [75l |76] . There are also several theoretical studies 
showing the hypothesized sensitivity of the 14 TeV LHC in the niA — tan (3 plane due to 
different search channels, for example Ref. m- In Figs. 1^ and we show the interplay 
between possible interesting signatures for H/A searches at the LHC and the GCE best fit 
regions in the — tan f3 plane, plotting contours of various branching ratios of interest for 
H/A searches at the LHC. To highlight the regions of interest, we overlay the values from 
Fig.|I] as dashed blue lines for Fermi spectrum (b) in Fig. |^and spectrum (d) in Fig. The 
gray shaded regions denote the current LHC exclusion limits from searches for H/A ^ r+r“ 
at the 8 TeV LHC {1-a and 2-a as labeled). In both figures, the two panels in the top rows 
show the branching ratios of the CP-odd Higgs: A —)■ (left) and A ^ Zh (right). The 

lower four panels display the branching ratios for H —)■ (middle left), H —)■ W^W~ 

(middle right), H ^ ZZ (lower left) and H ^ hh (lower right). 

The top row shows that both BR{A —)■ rr) and BR{A —)■ Zh) are a few percent through¬ 
out the parameter region of interest, with the former always comparable to or larger (in 
some cases, by more than an order of magnitude). We can understand this behavior by 
noting that due to the close to alignment conditions, the AZh coupling is very suppressed. 
Hence, despite the tan f3 enhancement of the gluon fusion production of A, we find that the 
rates for A ^ Zh are at least 2 orders of magnitudes smaller than the current exclusion 
limits [ZHl [79] and therefore unlikely to be probed even at the high luminosity LHC [75l [76] . 
Due to the absence of any other relevant decay modes, the decays to down-type fermions 
will still be the dominant decay modes and offer the best prospects for discovery of the 
pseudoscalar. 

For the heavier CP-even Higgs iL, in addition to the channel, there are non- 

negligible branching ratios into WW or hh despite being suppressed due to alignment (recall 
that, close to alignment, H ^ Hnsm)- These branching ratios are largest at low tan(3 below 
the top mass threshold, whereas Hr (FT —r+r”) is larger at higher tan/3. Note again that 
in the low tan (3 region, the main production of H is via gluon fusion, which is enhanced 
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FIG. 4: Predictions for LHC. Top-. BR{A r+ r") and BR{A ^ Z h). Middle: BR{H r+r") 


and BR{H W~^W~). Bottom: BR{H -A ZZ) and BR{H —)• hh). Dashed blue lines show 
values from fitting the GCE to spectrum-b, as seen in Fig. The colored contour regions (and 
bar on the right) are each plot’s respective bramiing ratio values. Shaded regions labelled 1-cr and 
2-f7 are the A/H —)■ r+T“ exclusion limits. 

























































FIG. 5: Same as Fig. but with values from fitting the Galactic Center excess to spectrum-d, 
as shown in Fig. 
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FIG. 6: Shaded contours denote values of 'R-ww Gray shaded regions bounded by solid black 
lines show 1-cr and 2-a exclusions by the H/A ^ searches from the 8 TeV LHC run (excluded 

above). The dashed blue lines correspond to contours of 


due to the large unsuppressed top coupling. We computed the rate of if —?• WW relative to 
the SM expectation, which is shown as colored contours in Fig. Current bounds 

on TZ-ww ^he level of 0.05 — 0.25 [HO], hence dedicated searches at the LHC could 

probe the GCE best-£t regions, particularly for uia ^ 350 GeV, where Tlww ^e within 
a factor of 10 of the current exclusion limit [751176]. 

For H/A heavier than about 350 GeV and low values of tan/3 (< 7), both the GP-odd 
and even Higgs bosons preferentially decay to top quark pairs. However, due to the large 
SM tt background, this is a very challenging signature for the LHG [8TI l82] : nevertheless, 
stronger sensitivity is expected at a 100 TeV collider [82]. The standard searches can 
probe regions with larger values of tan /3. 

It should be kept in mind that, in addition to these searches for heavier Higgs bosons, 
the good £t regions at low itia ^ 350 GeV also predict deviations in R^yy (see Eq. 
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for 


definition) at the 10% level or more, hence such deviations from SM-like properties of the 
125 GeV Higgs could be a stark signal of this scenario. All of the above search modes as well 
as the precision measurements of the 125 GeV Higgs are expected to improve substantially 
in sensitivity with the higher luminosity and energy of the 13 TeV LHG 
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B. Direct Detection 


Our predictions for spin-independent direct detection experiments are plotted in Fig.[^ 
which shows DM masses and spin-independent DM-nucleon (proton) direct detection cross 
sections compatible with the GCE (Fermi spectrum (b) in blue, spectrum (d) in red). We 
only show points with < 50 that are compatible with both the 2a A/H —8 TeV 


LHC constraints and 0.7 < Tlww — As discussed in Section IIIC, we see that DM via 
the pseudoscalar resonance corresponds to generic cross sections of C>(10“^^)pb, and these 
are comfortably safe from the existing XenonlOO [SI] and LUX jSS] bounds. A major fraction 
of the predicted parameter space can be probed with the next generation of direct detection 
experiments such as XenonlT and LZ [83] . We note that almost all points predicted from 
our £t lie above the neutrino floor and therefore a signal can in principle be detected. The 
green cross and star correspond to the best fit points from Fig.[^for spectrum (b) and (d) 
respectively. 


VI. SUMMARY 

To conclude, we summarize the main findings of this paper: 

• Recent reanalysis of GC background has found that the GGE could be consistent 
with annihilation of DM with much higher masses [ig [m dg 122]. This allows the 
GGE to be explained by the MSSM pseudoscalar resonance or “A-funnel”. We £t 
to two different dark matter spectra, Fermi spectrum (b) and (d) from [ig d2], and 
hnd that reasonable fits can be obtained while maintaining consistency with stringent 
constraints from collider searches, Higgs data, and direct and indirect detection. 

• For spectrum (b), the best £t region corresponds to 350 GeV < rriA ^ 450 GeV and 
tan/? < 6. This region can be probed with searches for H —)■ WW and ti resonance 
searches, uia ^ 250 GeV also gives reasonable fits but is incompatible with Higgs 
data. 

• For spectrum (d), there are two regions with reasonable fits to the GGE: 450 GeV 

^ ^ 600 GeV at tan/? < 8, and rriA ~ 300 and tan/? < 5.5. The former region 

can yield signals at the LHG in the A/H ^ tt oi ti resonance searches at the LHG. 
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FIG. 7: Dark matter masses, and spin-independent DM-nucleon (proton) direct detection 
cross sections, a^j, predicted by our fits to the Fermi GCE. Points compatible with Fermi spectrum 
b (d) are in blue (red); we have only plotted points with ^ 50 and compatible with collider 
and Higgs data (see text). The green cross and star correspond to the best fit points for spectrum 
(b) and (d) respectively. Fig. shows contour regions from fitting the galactic center excess to 
Fermi spectrum (b) and (d). Gurrent bounds (XenonlOO, LUX), the reach of upcoming detectors 
(XenonlT, LZ), and the neutrino background floor are also shown [83]. 

The latter region can also be probed with the same channels, and should also lead to 
measurements of deviations of the 125 GeV Higgs couplings from SM-like values. 

• The best £t regions for both spectra (b) and (d) predict spin-independent direct de¬ 
tection cross sections of C)(10“^^)pb for a 110 GeV < < 350 GeV neutralino. The 

entire region lies above the neutrino background, and the majority of the region is 
within reach of XenonlT and LZ (see Fig.[^. 

This exercise therefore leads to very sharp predictions for the next round of the LHG and 
direct detection experiments. Although the best fits obtained in this paper are noticeably 
worse than the best fit dark matter scenarios discussed elsewhere in literature, this highly 
predictive framework, coupled with the wide popularity of the MSSM, makes these results 
noteworthy. Even if the GGE turns out to be incompatible with the MSSM pseudoscalar 
resonance and is ultimately explained by some other (dark matter or astrophysical) phe- 
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nomenon, this study still serves as a valuable template for the interplay between existing 
collider and Higgs constraints and the indirect, direct, and collider signatures of the H-funnel 
region with a light pseudoscalar in the MSSM. 
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Appendix A: Parameters and Vacuum Metastability 
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Fig.|^ presents contour plots of the scanned parameters in the m^-tan/3 plane. The 
approximate check for vacuum metastability from Eq. 22 is shown in Fig. It is seen that 
the desired condition is satished (corresponding to the plotted ratio being less than 1) in 
most of the parameter space not ruled out by the 8 TeV LHC A/H ^ r+r“ bound. 


6 = |1-4m//m/| 
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FIG. 8: Contours of various input parameter values in the scan region. See text for details. 






























FIG. 9: Vacuum metastability requires this ratio to be approximately less than 1 gZ! , so we see 
that most of our points are compatible with vacuum metastability bounds. 
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